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We study a seesaw-type extension of the Standard Model in which the symmetry group is enlarged 
by a global (7(1). We introduce adequate scalar and fermion representations which naturally explain 
the smallness of neutrino masses. With the addition of a viable scalar Dark Matter candidate, an 
original scenario of leptogenesis emerges. We solve the relevant set of Boltzmann equations and 
show how leptogenesis can be successfully implemented at the TeV scale. The constraints on the 
scalar mass spectrum are derived and the Dark Matter phenomenology is discussed. 
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I. INTRODUCTION 



Now that we entered in the LHC era, the Standard Model (SM) of elementary particles can be definitively tested. 
Until now, the SM has been extremely successful, as no strong signals of new physics have been observed so far 
at particle accelerators. However other experiments have long-time evidences for the need of extensions of the SM 
particle content. Neutrino oscillations are the prime among them on the particle side, but the compelling gravitational 
evidences for the existence of Dark Matter (DM), as well as the observation of a matter-antimatter asymmetry in the 
Universe all call for new physics. 

From neutrino oscillation experiments we know that at least two neutrinos should be massive with an overall mass 
scale m v constrained by different observations: m v < 1 eV. More precisely, experiments with solar, atmospheric, 
reactor and accelerator neutrinos [I]-[1Q] set two mass scales in the theory, Atoq and Am}, which drive the solar and 
atmospheric neutrino oscillations, respectively [llj : 

Am| = (7.59±gJS) x 1CT 5 eV 2 , Am\ = (2.43 ±0.13) x 10~ 3 eV 2 . (1.1) 

Moreover, these experiments show that flavor neutrino mixing, described in terms of the PMNS [T2T[T^] matrix, is 
characterized by two large mixing angles, 6*12 and 6*23, and a small one, #13 [T5] . 

On the cosmological side, the matter content of the Universe has been measured with precision by WMAP [16] . 
The resulting Dark Matter and baryon number densities, JJdm and flB, are 

fi DM = 0.229 ±0.015, fl B = 0.0458 ±0.0016. (1.2) 

Several gravitational observations confirm the existence of non-baryonic matter [17] , which is not accounted for in 
the SM. New physics extensions are then necessary and various viable DM candidates exist [17]. However, the real 
nature of DM is still elusive, as no direct proof has been observed - or firmly confirmed - so far [18] -[21]. The 
measurement by WMAP of the baryonic matter content of the Universe is in agreement with the value predicted by 
Big-Bang Nucleosynthesis from the observations of the primordial abundances [22 . However, an excess of baryons 
over antibaryons is observed, and the standard cosmological scenario fails to explain this Baryon Asymmetry of the 
Universe (BAU). Particle physics extensions of the SM are advocated to justify this: in relation with neutrino masses, 
the leptogenesis scenario [23 [21] constitutes one of the most elegant solutions. 

In this paper we study a minimal extension of the SM in which it is possible to address, in a consistent way, the 
three puzzles listed above. The model is based on a global U(1) B _? symmetry, which is spontaneously broken below 
the electroweak symmetry breaking (EWSB) scale. The L charge is a generalization of the usual lepton number L, 
a,s L — L for the SM particles. The light neutrino masses are explained within a seesaw framework [25], through the 
introduction of a SM singlet Dirac fermion Nj), together with three Brout-Englert-Higgs scalar particles: two SU(2)w 
doublets Hi t 2 and a SM singlet H 3 , which drive the EWSB by acquiring non-zero vacuum expectation values (vevs). 
All these extra degrees of freedom are charged under the global U(l) B _i symmetry. In e.g. [26], neutrino masses 
were generated in models with similar scalar spectrum and/or based on a (spontaneously broken) global symmetry, 
although in different physical frameworks. In our scenario, when the seesaw scale is set in the TeV-range, such a 
particle content provides a UV-completion of the inverse-seesaw mechanism of neutrino mass generation [27] . 

Nevertheless, with just this particle content, neither the observed amount of baryon asymmetry nor the Dark Matter 



abundance, eq. (1.2), can be accounted for. 

In order to solve also these two important issues, we complete the model by introducing a Majorana neutrino A3 
and a complex scalar S. Both particles are SM singlets, although S is charged under the global [7(1) B _£. The particle 
content of the model is summarized in Tab. |TJ together with the U(l) B _i quantum numbers of the fields. The new 
scalar S provides, after the breaking of U(1) B _£ ) a natural Dark Matter candidate, whose stability is guaranteed by 
a remnant Z% symmetry. 

It is remarkable that the introduction of S allows a TeV scale scenario of leptogenesis. Indeed, as the Majorana field 
A3 couples to Ne> and S, the out-of-equilibrium CP-violating decays of N 3 can generate a number density asymmetry 
in Njj and 5, resembling the standard thermal leptogenesis mechanism in the type I seesaw extension of the SM. 
However, in the present case leptogenesis is implemented in two steps: first an asymmetry in Np and S is generated 
by the decays of A3; in a second phase, the Dirac neutrino asymmetry is transferred to SM leptons by sufficiently 
fast neutrino Yukawa interactions. The latter set a link between successful leptogenesis and viable neutrino mass 
generation via the seesaw mechanism. Finally, as in standard leptogenesis, non-perturbative sphaleron effects partly 
convert this lepton asymmetry into a net baryon number [2"B] . 

In Section [IT] we discuss neutrino mass generation through the (inverse) seesaw mechanism. In the subsequent 
Section [Hi] we tackle the problem of the BAU and study the constraints on the parameter-space of the model imposed 



by successful leptogenesis. The computation of the CP asymmetry and the set of coupled Boltzmann equations 



governing the number density evolutions are reported in the final appendices. In Section IV we discuss the scalar 
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TABLE I: Charge assignment of the fields. 



sector of the theory, deriving the mass spectrum and corresponding constraints. In Section|V]we study the possibility 
of having a viable Dark Matter in the model and comment on the possible observation of DM in direct detection 
experiments. Finally, in the last section we summarize the main results of the paper. 



II. NEUTRINO MASSES WITH A GLOBAL U(1) B _ L 

An effective Majorana neutrino mass term is generated below the EWSB scale from the following part of the 
interaction Lagrangian: 

- £ int D MN D N D + (y\ N D H{ l t + y{ N D H\ t d + ^= H 3 N D N C D + h.c)j , (2.1) 

where £, = (viL,£iL) T (i = e, /i,r), N D = CN D and Hk = —iaiHT (k — 1,2). 1 The coupling constant a and the 
neutrino Yukawa couplings y\ 2 are complex parameters. As we will see in the following, the phase of a plays a crucial 
role in the generation of the CP asymmetry necessary for the production of the observed amount of BAU. 



The terms reported in the Lagrangian (2.1 1 provide a dynamical realization of the inverse seesaw mechanism |27j 
for the generation of neutrino masses in the case the mass of the Dirac field Nd is taken in the TeV-range. More 
specifically, in our scenario the standard lepton charge L is explicitly violated by the interactions involving the 
couplings y l 2 and a. Consequently, we expect that the active neutrino masses, generated through the (inverse) seesaw 
mechanism, do directly depend on these parameters. The model, in this minimal form, predicts two massive and one 
massless active neutrinos. 

The seesaw mass scale M is a free parameter of the theory and can assume arbitrarily large values above the EWSB 
scale. However, in the following we will be mostly interested in the case where M is taken at the TeV scale. At 
energies much smaller than M, Nd is integrated out and we get at second order in 1/M the (B — L)-conserving 
effective Lagrangian: 2 



a. ^ - (?£;) {ale,) + 4fg p; si) (m) * 



3 

(2.2) 



(fiH*) (hU % ) H 3 + h.c, 



V2M 2 

where the sum over the flavor indices i and j is understood. When the neutral components of the scalar fields Hk 



(k = 1,2,3) take a non-zero vev, the operators in (2.2 1 generate a Maj oran a mass term for the flavor neutrino fields 

T 

Vi l . Indeed, taking (Hi) — (0, Vi/y/2) (i = 1,2) and (^3) = vs/y/2 in (2.2|, we obtain the neutrino mass Lagrangian 



C m „ = -^n.n% u v-l + h.c. , (2.3) 

where j/ l = (^eL^tiL^rL), = CviJ and 

1 \ ( i 3 , i j i j *ViV 3 i j v 2 v 3 \ viv 2 ,„ ,x 

M l3 = (Vi Vi + Vl Vl -Vlvia —-Ma—)—. (2.4) 



1 C is the usual charge conjugation matrix of Dirac spinors. 

2 We do not include flavor kinetic mixing terms in the Lagrangian ( 2 . 2 i , which arise by dimension 6 effective fermion operators. 
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The masses of the two active neutrinos are given by 
1 



v 2 2 v l 2 * ^ v 2 f v 2 9 v l 9 ^ w 2 \ . o Wo » 



,(2.5) 



where we define y 12 = y\y e 2 + y{y^ + y[y T 2 , y k = \/(y e k ) 2 + (Vk) 2 + (vl) 2 ( k = 1 , 2 ) and 7712 = 
V (2/12/2 ~ y^Vi) 2 + (2/12/2 — ytVi) 2 + (2/12/2 — 2/2 2/i) 2 - As usual in two-Higgs doublet models, the vevs of the two 
scalar doublets, V\ and «2, are related to the EWSB scale: y/vf + «| = i> — 246 GeV. As explained in Section |iv| 
the hierarchy among the Higgs vevs is tightly constrained in our model, in particular from the presence of a massless 
Goldstone boson associated with the spontaneous breaking of the global U(l) R _f. phenomenological constraints en- 
force v 2 "C t>i,3, and by convention we impose v 3 < v. As we will see in Section [IV[ this hierarchical pattern is easily 
realized in the model. Typically, for \a\ « 0.01, Mwl TeV and a scalar spectrum with v 2 ~ 10 MeV, U3 ~ 100 GeV, 
the neutrino Yukawa couplings are 1 2/1,2 1 ~ 10~ 4 . 

The Yukawa interaction aH 3 No Np generates after EWSB a small Majorana mass term for the two chiral com- 
ponents of the Dirac field Njj, which is then split into two quasi-degenerate Majorana fermions: they behave as a 
pseudo-Dirac pair |29j-|31j. with a mass difference of the order 2u 3 |a|. Such scenarios have been studied in detail in 
[52] . where it was shown that a high-level of degeneracy prevents the Majorana nature of these states to be observed 
at colliders, LHC included. Indirect signals of TeV scale pseudo-Dirac neutrinos coupled to charged leptons can in 
principle be observed both in lepton flavor violating processes, e.g. charged lepton radiative decays t\ — > lj 7 and 
fx — e conversion in nuclei, and in experiments searching for lepton number violation, such as neutrinoless double 
beta decay processes. For these processes, the contribution of the heavy neutrinos to the decay rate may be rel- 
evant/dominant in the case of M « (100 — 1000) GeV, \a\v 3 /M rj 10 -3 — 10~ 2 and for sizable neutrino Yukawa 
couplings, 1 2/1,2 1 ~ [33]. 

Finally, we remark that the coupling a is not strictly required in order to obtain two massive neutrinos, whereas 
the introduction of j/2 is mandatory. Actually, one can show that y\ and y 2 are also sufficient to fully reconstruct the 



low-energy neutrino data, up to a normalization factor [34 j . From eq. (2.5) we get the following relation 



|r?i 2 | v 1 v 2 ^ = 2 (Am| Am|) 1/4 • (2-6) 
This equation clearly shows that for y 2 — or for (y 2l y%, y 2 ) aligned with (yf , y^,y\), \rj\ 2 \ — and only one neutrino 



is massive, in contradiction with neutrino oscillation data. Barring accidental cancellations, eq. (2.6| implies 



„ s ( M \ /10 MeV\ 

III. TWO-STEP LEPTOGENESIS 

Before discussing how the baryon asymmetry is generated in our scenario, let us briefly recall the standard picture 
of leptogenesis, based on the type I seesaw extension of the Standard Model. For a detailed discussion, see [3S] and ref- 
erences therein. In the standard scenario, at least 2 massive right-handed (RH) neutrinos, which are SU(2)w x U(l)y 
singlets, are introduced and couple to lepton doublets through Yukawa interactions. These singlets are Majorana 
fermions whose mass Mr is not related to the electroweak scale and can assume arbitrarily large values. The RH 
neutrinos evolve together with the SM particles in a hot but expanding Universe; when the temperature drops down 
below Mr, they start to decouple and decay out-of-equilibrium in both leptons and antileptons. If CP is violated in 
these processes, a non-zero asymmetry is produced, which is subsequently converted into a net baryon number by fast 
sphaleron interactions. The latter are non-perturbative effects, in thermal equilibrium above the EWSB scale up to 
temperatures T < 10 12 -10 13 GeV [35]. Several interactions should be considered for an accurate determination of the 
efficiency of leptogenesis in producing a baryon asymmetry. Spectator processes play an important role in modifying 
the production/depletion mechanisms, most notably by spreading the lepton asymmetry into different species. 

In the present case, given the particle content and the charge assignment listed in Tab[lJ the interaction Lagrangian 



receives, besides the operators of the seesaw sector in eq. (2.1), contributions from the extra Majorana field A3 and 
the scalar S: 

1 / ,/< 

A/f J\T ATC 1 I „ C* AT AT ^ o2 i 



-£ int D Ms^*^ + N * N Z + [gSN D N 3 - ^S 2 H* 3 + h.c. ) , (3.1) 
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where M3, fx" and g can be set real by a redefinition of the phases of N3, S and H3. We impose A3 to be heavier 
than No and S. 

In this model, the generation of a baryon asymmetry proceeds in two different phases. In a first phase, which is 
similar to the standard leptogenesis scenario, an asymmetry in Dirac neutrinos Njj and in S is generated by the 
out-of-equilibrium decays of the Majorana field A3. As we describe below, the CP asymmetry in N3 decays is only 
possible after the introduction of S, carrying the same B — L quantum number as Njj. 

Besides decays and inverse decays, several scatterings affect A^ and S asymmetries. All these interactions conserve 
the total B — L charge. In a second phase, owing to the neutrino Yukawa couplings, the produced Nd and S 
asymmetries are transferred and reprocessed into a lepton asymmetry. In this second phase, the sphaleron processes 
partly convert the so produced lepton asymmetry into a final baryon number, as in the standard picture. 

This model can thus be viewed as the SM augmented with a second Higgs doublet, combined with a hidden sector 
composed of the fields N3, S and H3. The two sectors share a conserved B — L charge through the Dirac neutrino Nr>. 
In that extent, the role of the neutrino Yukawa couplings is central both in the generation of light neutrino masses 
and in the production of a BAU, in agreement with observations. 



A. The CP asymmetry ecp 



In the standard leptogenesis scenario a CP asymmetry is generated by the interference between the tree-level and 
the one-loop corrections to the decay amplitude of the heavy Majorana neutrinos (2~il |3"T] . owing to the presence of 
at least two heavy states. In our case, with only one heavy neutrino Ad, no CP violation is produced in Njj decays. 
On the other hand, a non-zero CP asymmetry can be generated by the addition of A3 and S, from the interference 
between the tree-level and one-loop correction to N3 decay amplitude, whose Feynman diagrams are depicted in Fig. [I] 

The detailed computation of the CP asymmetry in N3 decays is provided in Appendix [Bj We report below the 
resulting expression in the limit M 3 ^> M, fig- 



Im(a) n" 
16tt W3 



(3.2) 



Despite of the fact that A3 decays depend on the coupling constant g, the latter being a real parameter does not 
enter in the expression of ccp, cf. eq. (B13|. The only source of CP violation relevant for leptogenesis is the phase 
of the complex parameter a in the Lagrangian (3.1). It is remarkable that, in contrast to the standard leptogenesis 



scenario, there is no direct dependence of ecp on the neutrino Yukawa couplings 2/1,2- Still, a connection between the 
leptogenesis CP- violating phase and the light neutrino masses exists and is actually provided by the imaginary part 
of a. We remark that the parameter /1" in ( |3.2[ ) enters in the mass splitting between the real and imaginary parts of 
S (cf. eq. (5.1 )) and therefore determines which is the DM candidate of the model, as shown in Section [V] Provided 
fj," is not too much suppressed compared to M3 and the phase of a is different from zero, ecp takes sizable values. 
We typically have: 



ecp 



-2 x 10" 



6 (j^L\ f 10TeV ^ 

lGeVy \ M 3 J 



Im(a) 



(3.3) 



6 



B. Asymmetry productions 

We discuss now the salient aspects of leptogenesis in our scenario. We eventually distinguish between two stages 
of production, but we shall emphasize that these stages are not necessarily consecutive and may occur in the same 
temperature range. 



N D Nn _ N D N D _ _ N D N D 




FIG. 2: Feynman diagrams of the ANd = AS = 2 scatterings. 




FIG. 3: Feynman diagrams of the ANd = AS = 1 scatterings. 



x S H ly x S H lf 
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FIG. 4: Feynman diagram of the S self-annihilation. 
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g = 1CT 5 , |o| = 1 g = 1CT 3 , |tt| = 1 g = IO -3 , |a| = 1CT 2 




icr 2 icr 1 i io 1 io 2 icr 2 icr 1 i io 1 ia 2 icr 2 icr 1 l io 1 io 2 



FIG. 5: Processes relevant in the first stage of leptogenesis: thermal density rates as function of z = Mz/T, for M3 = 50 TeV, 
M = 10 TeV and fj." = 1(100) GeV, plain (dashed) curves. In black is reported the total decay rate of N3. The purple curves 
stand for the sum of the (non-resonant part of the) diagram a) and diagram b) in Fig. [2] The blue, orange and red curves 
correspond, respectively, to the processes a), b) and c) shown in Fig. [3] 



1. First Stage: processes at 0(a 2 ), 0(g 2 ), 0(g 2 a 2 ) and 0(g 4 ). 

We list below the processes relevant in the first step, where the asymmetries in S and are created. 3 Further 
details are given in Appendix [Cj 

• Decays and inverse decays of N3: N3 — > Nrj S , N d S (see Fig. [lj. 

• AN D = AS = 2 scatterings: N D Sd- N D S and N D N D o S S (see Fig. [2). 

• AN D = AS = 1 scatterings: N D N 3 o H 3 S, N D H 3 o N 3 S and N D S <-» iV 3 H 3 (see Fig. |3j|. 

• S self-annihilation: S S o Hi H 2 (see Fig. U 



Notice that the last process depends on interaction terms reported in the scalar potential of the model (see eq. (4.3)). 
However, it turns out to be numerically irrelevant, so we disregard the effect of this term in the following. 

We display in Fig. [5] the interaction rates j eq of some of the up-listed processes as function of the parameter 
z = M 3 /T, where T is the temperature of the plasma. These rates are normalized by H(z)n e ^ (z), except for the 
ANd = 2 rates which are normalized by H(z) n e ^ D (z), as they only act as damping terms. For illustration we fix 
M 3 = 50 TeV, M = 10 TeV and we choose representative values of g and |a| for the different panels. 5 

We represent in each plot by straight (dashed) lines the computed rates assuming //' = 1 (100) GeV. For //' = 1 GeV, 
the cross-sections of the ANd — 2 scatterings a) are dominated by their s- and u-channels and scale as 0(g 4 ). The 



3 We denote by AX the absolute variation of the X particle number density. 

4 In Fig. [5] only the off-shell part of the AJVp = 2 diagrams a) is shown, as its on-shell part equals 715/4, -yjj being the total decay rate 
of N3. 

5 For definiteness, in the following numerical evaluations we set the phase of a to its maximum value a = — i [a|. 



n" = 1 GeV 



It" = 100 GeV 



\a\ 10 




FIG. 6: Efficiency rji of production of Nd asymmetry in the first stage of leptogenesis, cf. eq. ( |3.4[ ), as a function of g and \a\ 
for Ms = 50TeV, M = lOTeV and /i" = 1(100) GeV in the left (right) panel. 



ANd = 2 scatterings b) on the other hand are governed by their s-channcl and are proportional to \a\ fjf . The 
ANo = 1 processes a), b) and c) of Fig. [3] are dominated by their respective t-, s- and it-channels and therefore scale 
as 0(g 2 \a\ 2 ). For larger values of //', e.g. fi" — 100 GeV, the ANo — 1 processes b) and c) get sizable contributions 
from their i-channels (cx g 2 /x" 2 ) which dominate over the other channels for small values of a, as can be seen in the 
right panel of Fig. [5] The different interaction rates where evaluated using the packages FeynArts [38] and FormCalc 
[39] . To this end, we implemented our model, £ int eq. (Al), via FeynRules [40] . 



The various interactions considered above control the amount of No and S asymmetries produced during the first 
stage of leptogenesis. As the lepton asymmetry -and finally the baryon asymmetry- mostly depends on the amount 
of No asymmetry produced in the first step, it is useful to introduce an efficiency factor r/i defined through: 



YaNo (*tr) = £cp m YZ (T » M 3 ) . 



(3.4) 



In this parametrization, Yx indicates the comoving number density of X, while z tr ~ M 3 /M approximately marks 
the transition between the first and second stage: for z > Zt r , i.e. T < M, No decouples from the plasma and decays 
into leptons and antileptons. 

Given the numerous interactions considered above, the derivation of an analytic expression for the efficiency factor 
rji is quite challenging. Nevertheless, we perform a numerical evaluation of r]i by solving the set of Boltzmann 
equations reported in Appendix [Cj The resulting efficiency is shown in Fig. |6j where iso-contours of 771 in the g — \a\ 
plane are displayed, for M 3 = 50 TeV, M = 10 TeV and p!' = 1 GeV (100 GeV) in the left (right) panel. 

We first consider the case of small fi", left panel of Fig. |6j In this case, the AA^d = 2 scatterings are typically 
smaller than the decays and inverse decays, as shown in Fig. [5] Depending on the value of a, the ANo = 1 scattering 
rates 7^ (k = a, b, c), may be in equilibrium when the No asymmetry is produced. This occurs if 



7^ 



H(M 3 ) 



> 1 



a x 



Go- 6 ) 



> 1 



(3.5) 



As can be seen in the left panel of Fig. [6] the efficiency rji strongly depends on whether the ANo scatterings are in 
equilibrium at T ~ M3 or not. In the case their rates are not fast enough, i.e. if the condition (3.5) is not satisfied, 
the production of No and S asymmetries is mostly driven by decays and inverse decays of A 3 . This situation is very 
similar to the standard leptogenesis scenario, when AL — 1 scatterings are neglected. Therefore, we expect that larger 
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values of the coupling g increase the washout effects. The strength of N3 decays and inverse decays can be expressed 
in terms of the washout parameter Kd 

v - T Nz „ ( 9 \ 2 {50TeV\ 



For g > 10 6 and M3 ~ O(10) TeV, decays and inverse decays act in a strong washout regime, where the efficiency is 
approximately given by [41 : 

0-4 

* ~ K D log(^) ' (3J) 

For smaller values of decays and inverse decays act in a weak washout regime, and the efficiency scales as K D [41j . 
in the case where the abundance of N 3 is vanishing at high temperatures. 



In the opposite regime, when the ANd = 1 scatterings are fast enough and the condition eq. (3.5) is satisfied, an 
initial (anti-)asymmetry is produced at earlier times, due to the CP violation in scatterings, which is discussed in 
Appendix [c| From Fig.[6j we can distinguish two relevant cases, according to the values of a and g. For \a\ ~ 10~ 6 /(/, 
the ANd = 1 scatterings essentially act as source terms, producing Nd and S asymmetries, thereby increasing the 
efficiency r\\. This effect is manifest in the diagonal of the left plot of Fig. [6} Conversely, for larger values of a, the 
ANd — 1 scatterings act as damping terms and increase the washout of the asymmetries. The resulting efficiency is 
therefore highly reduced. 

The case of a larger fj," is depicted in the right panel of Fig. |6j where we fix fx" = 100 GeV, while M 3 and M 
assume the same values as before. As already stated, in this case the ANd = 1 scatterings b) and c) pick-up sizable 
contributions from their corresponding i-channels, and are enhanced for relatively small values of a, compared to 
the fj," = 1 GeV case, as can be seen in the right panel of Fig. [5] For |a| < 0.1, the efficiency depends essentially 
on g, since the scatterings are negligible with respect to the decays and inverse decays, and then r/i behaves as in 



eq. (3.7). However, for small values of g, g < fewlO -5 , ANd = 1 scatterings become competitive with decays and 
inverse decays, both in the generation and in the washout of Nd and S asymmetries. For values of a of order one 
and small g, the scatterings mainly act as source terms in analogy with the ji — 1 GeV regime, thus increasing r]i. 

We see that in this first stage, the efficiency of Nd and S asymmetry production can be close to its maximum 
possible value in a large region of the parameter-space. However, this does not guarantee a successful leptogenesis, as 
this asymmetry should be transferred efficiently to leptons. 



2. Second Stage: processes at 0(2/12), 0(g 2 2/1,2)- 



We now concentrate on the second step of leptogenesis: the transfer of Nd asymmetry to the lepton doublets. Once 
a lepton asymmetry is generated, the sphaleron processes which are active at the leptogenesis epoch convert part of 
it into a non-zero baryon number density. The second stage ends at the freeze-out of the sphalerons, that may occur 
before or right after EWSB [42]. 

We report below the main A£ = 1 processes which participate in the lepton charge transfer mechanism: 

• Decays of Nd, which are either L-conserving, Nd IHi, or L- violating, Nd —> £ H 2 . 

• Scatterings on top quarks: the s-channel Nd~£ *h> tq^ and the i-channcls NdQ3 (t) O lt(cfe). These processes 
are mediated by the exchange of the Higgs doublet Hi and correspond to the AL = 1 scatterings in standard 
leptogenesis. Notice, however, that in our case lepton number is conserved. 

• Scatterings on iV^: N3 S ■(-> £ Hi and N3 S -H- £H 2 which are mediated by Nd- A CP asymmetry emerges from 
these processes, as shown in Appendix |C| 

We do not include the scatterings involving gauge bosons in our evaluation of the baryon asymmetry. However, we 
do not expect these processes to have a quantitative impact. Indeed, they cannot act as a source term for the lepton 
asymmetry since no CP violation is possible in this case, in contrast to the standard leptogenesis scenario |43j . In 
addition, they tend to equilibrate the lepton and Nd number densities, like the scatterings on top quarks considered 
above. Actually, it is shown in references [H]-[1S] that these processes have comparable rates. 

The lepton doublet can also participate in AL = 2 A^-mediated scatterings, similarly to the standard leptogenesis 
case: £ Hi f-> £Hi and W f> Hi H 2 - In this case the scattering rate is proportional to both the neutrino Yukawa 
couplings, j/i and 2/2 • In a democratic scenario, that is for \yi\ w |y 2 |, provided the constraints from active neutrino 



masses, eqs (2.5) and (2.6) are satisfied, such AL = 2 scatterings are usually in equilibrium at the leptogenesis time. 
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g=icr 3 , oil = \H = 10^" 




FIG. 7: Processes relevant in the second stage of leptogenesis: thermal density rates y eq as function of z = M3/T, for 
M3 = 50 TeV, M = 10 TeV. The blue curves stand for decays of No- The orange curves stand for the (sum of s- and f-channel) 
At — 1 scatterings on top quarks, while the purple ones represent Al = 1 scatterings on Y3. The green line corresponds to the 
total AL — 2 scattering rate. 



They are however greatly suppressed compared to the A£ = 1 scatterings and turn out to be numerically irrelevant, 
as illustrated below. 

In this second stage, all interactions depend on the neutrino Yukawa couplings y\ and yi- In the limit where these 
couplings are zero, no lepton (doublet) asymmetry can be generated as basically both No and S decouple from the 
SM sector. This clearly implies a lower bound on the values of y\ and 2/2 ■ 

Let us discuss this bound, independently of the constraints from low-energy neutrino masses, as it sheds light on 
how this second stage works. To this end, we represent in Fig. [7] the processes relevant in the second step for the same 
set of parameters used in Fig. [5} M 3 = 50TeV, M — lOTeV and g = 10~ 3 , and for \yi\ = |y 2 | = 10 -4 . We represent 
in Fig. [7]by plain (dashed) curves the rates normalized by H(M), where X = No (I), acting as source (damping) 
terms. The blue curves correspond to No decays, both L-conserving and violating as \y\ \ = \yi\- The orange (purple) 
curves are related to scatterings on top (N3), while the green line stands for the AL = 2 processes. Scatterings on 
N 3 and AL — 2 interactions are clearly sub-dominant and can be neglected. 

A lower bound on y\ can be derived by demanding that the scattering rates on top quarks, denoted by r y% D , are in 
equilibrium at T <~ M, when acting as a source term for the lepton asymmetry: 



Therefore, provided yi is large enough, the L-conserving scatterings are in equilibrium and can transfer the Njj 
asymmetry to the lepton doublets. A similar lower bound arises for yi from the corresponding AL = 1 scatterings 
with gauge bosons. 

The main source of lepton asymmetry production may originate just from the decays of No- Let us suppose, 
indeed, that the lower bound on yi, eq. (3.8), is not satisfied. Still, as we see from Fig. [7j decays dominate over the 



scatterings at T ~ M. For these decays to be effective in redistributing the No asymmetry to lcptons, the Dirac 
neutrino should be heavy enough, say M > 10 x T sp h, and the following condition should be satisfied: 



T Nd >H(M) | yii2 |>6xl0- 7 y i5 ^. (3.9) 

In summary, for neutrino Yukawa couplings smaller than the bound above, the lepton number asymmetry production 
is not efficient, if only condition ( |3.9[ ) is satisfied, almost all No decays to leptons, and we expect that at the end of 
the second stage, the lepton asymmetry equals the amount of No asymmetry produced in the first stage. For larger 
Yukawa couplings satisfying eq. ( |3.8| ), Y&i{z tr ) Y&]y D (z tr ) at the end of the first stage, so at the end of the second 
stage, Y Ae (z sph ) w 2Y AND (z tr ). 
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FIG. 8: Influence of scatterings on the transfer of Nd asymmetry to a lepton doublet asymmetry. See the text for details. 



The case of a light Dime Neutrino 



An interesting case is realized when the Dirac neutrinos are so light that they don't have enough time to decay before 
the freeze-out of the sphalerons. Deman ding that the scatterings with quarks are in equilibrium, at least slightly before 



the sphalerons decouple, the condition (3.8) is changed to \yi\ > 10 -6 . Therefore y\ should be at least larger than the 
electron Yukawa coupling. We illustrate this remarkable case in Fig. [5J where we fix M = 200 GeV and we consider 
two sets of values for \yi\ and 1 2/2 1 : i) |t/i| = j 3/2 1 = 5 X 10~ 4 (red curves) and ii) \yi\ = I2/2I = 5 x 10~ 7 (black curves) 



Notice that the latter case may hardly be compatible with constraints from neutrino masses, eqs ( |2.5[ )-(2.7l. The 
dashed and plain curves correspond to the lepton doublet and Nn asymmetries, respectively, and all the asymmetries 
have been normalized to Yn d (ztr), where z tr ~ M3/M is indicated by a blue band. We impose a sphaleron freeze-out 
at around T sp h ~ 130 GeV, which is represented by the gray band. 

We see from Fig. [8] that while Njj asymmetry is almost unaffected by the Yukawa hierarchy, in case ii) leptons do not 
equilibrate with Nd as scatterings are out-of-equilibrium, while in case i) Yai ~ Yan d at temperatures well above 
M. 

Provided that the neutrino Yukawa couplings are sufficiently large, an asymmetry in Nd will be always transmitted 
to the lepton sector, regardless of the Dirac neutrino mass: we can therefore asset that no lower bound can be derived 
on M from leptogenesis. 

In conclusion, once light neutrino mass constraints are applied, a lepton asymmetry is efficiently produced. A 
successful leptogenesis then only relies upon the amount of Nd asymmetry produced in the first stage. 



C. Successful leptogenesis 



In the former subsections, we analyzed the conditions under which a Nd asymmetry is efficiently produced during 
the first step of leptogenesis, and subsequently transmitted to the lepton doublets. Through the sphaleron processes, 
this lepton asymmetry is partly converted into a baryon number density. The sphalerons violate both lepton and 
baryon numbers, but conserve B — L: it is therefore more convenient to evaluate the B — L asymmetry. Given the 
different processes in thermal equilibrium during leptogenesis era, the final baryon asymmetry reads: 

Yab = y YA(B-L){z sp h) ■ (3.10) 



The derivation of eq. (3.101 is given in Appendix [d] We illustrate in Fig. [9] the evolution of Nd, S and baryon 
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If = 1 GeV ji" = 100 GeV 




FIG. 9: Evolution of number density asymmetries in function of z. M s = 50 TeV, M — 10 TeV, \a\ = 0.1 and g — 10 5 
are fixed. The black plain (dashed) curves represent ecp x Y^?* ■ In dotted-blue, dot-dashed orange and plain red are shown 
Yan d , Yas and Yab respectively. 
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FIG. 10: Successful leptogenesis. Region of the parameter-space in the |a|-<? plane providing a final baryon asymmetry (not) 
compatible observations, (yellow) black points. We fix M 3 = 50 TeV, M = 10 TeV, fx s = 100 GeV and fi" = 100 GeV. 



asymmetries against z for typical values of the parameters: M 3 — 50 TeV, M = 10 TeV, /is = 100 GeV and fj," = 1 
or 100 GeV, left or right panel respectively, for fixed values \a\ = 0.1 and g = 10 -5 . For such values of \a\ and g, 
we see from Fig. [6] that the efficiencies rji are quite similar, that is ?/i ~ 0.1. However, for p," = 100 GeV the CP 



asymmetry is ecp — 4 x 10 , two orders of magnitude larger than for pc = 1 GeV, cf. eq. (3.2 1, and so the baryon 
asymmetry in the former case will be bigger. Indeed, for fj," = 1 GeV, Fab ~ 1-6 x 10~ n while for //' = 100 GeV, 
Yab ~ 10~ 9 . These values should be compared with the measurement of WMAP [IB]: 

= (8J7 ±o.21) x 10~ n . (3.11) 
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In Fig. 10 we made a scan over the two parameters a and g: the black points represent values of Yab compatible with 
observations. 6 As we see, a successful leptogenesis is easily realized in our scenario, provided the CP asymmetry is 
big enough, that is € CP > 3 x 10~ 7 , and the washout processes do not suppress the N D asymmetry in the first stage, 



IV. THE SCALAR SECTOR 



Given the charge assignment of the scalar fields in Tab. [TJ the most general scalar potential Vsc invariant under 
SU{2)w x U(1)y x [£^(1)b_£,] can be written in the following form 



Vsc = Vsb + Vdm , 

where Vsb and Vdm denote the symmetry breaking and dark matter scalar potentials, respectively: 



(4.1) 



V SB = 



-Mi Hi H x + \ x (Hi H,) 2 - fi 2 2 Hi H 2 + A 2 {Hi H 2 f - H* 3 H 3 + A 3 (H*H 3 ) 2 
/sia Hi H X H\ H 2 + k' 12 H\ H 2 Hl H x + /ti 3 h\ H x H^H 3 + n 23 h\ H 2 H* 3 H 3 

^= (hIh^ + hI^h; 



V 



DM 



/4 s* s + a s (s*sy + j"x h\ HiS*s + t 2 h\ h 2 s*s + j" 3 h;h 3 s*s 



h s 2 hI h 2 + h* s* 2 h\ h x - ^(s' 2 h; 



S* 2 H 3 



(4.2) 



(4.3) 



Through rotations of the scalar fields, all parameters but h can be made real, while the dimensional parameters are 
assumed positive. The parameter h is in general complex, but we will assume in the following that h is real. 

The two scalar doublets Hi 2 and the complex scalar singlet H 3 are responsible for the breaking of SU(2)w x 
U(l)y x down to U(l) em x \Z 2 \. Given the charges of H 2 and H 3 , the discrete Z 2 emerges as a remnant 

symmetry of the global U(l) B _r after EWSB. Among the ten real scalar degrees of freedom, three of them are eaten 
through the Higgs mechanism, leaving a spectrum of seven physical scalars: two charged particles, H ± , two CP odd 
neutral scalars, A and the massless Majoron J ^45], and three CP even neutral scalars, h°,H° and Ha- We derive 
in the following subsections some constraints on the scalar sector parameter-space. An exhaustive phenomenological 
study, although of great interest, is beyond the scope of this work. 

The minimization of the scalar potential with respect to Hi, H 2 and H 3 vevs enforces three tree-level relations, 
that we use to define the quadratic terms Indeed, by parametrizing the Brout-Englert-Higgs fields and S as 



with 



we get the extremum conditions 



dVsB 
dv. 



Hk 
H 3 



= H 



Vk + h k +i a k 



v 3 + h 3 + ia 3 



V2 



, fc = l,2, 
Sq + i Si 



(Hi) = ~ and (S)=0, 
V2 



VjVk}J_ 

2vi 



i- J? k — Ij 2, 3 , 



(4.4) 
(4.5) 

(4.6) 
(4.7) 



where ki 2 = k\ 2 



' 12 



and Ru = 



elsewhere. The extremum obtained in (4.7) is an absolute minimum provided the 



Hessian of Vsb is positive definite. Boundedness from below of the scalar potential requires the quartic couplings Afc 
to be positive, as well as a non-trivial relation among the couplings. Notice that, since both Hi and H 2 are charged 
under SU(2)w x U(1)y, they both contribute to the masses of the SM gauge bosons. 



6 Actually, for the sake of illustration, we enlarge the required range, demanding 3 X 1CT 11 < Y AB < 3 X 10 
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Among the numerous parameters of Vsc> it is worth to emphasize the role of the trilinear coupling \j! '. In [4T1 148] . 
a two-Higgs doublet model was built invariant under a U(l) global symmetry, explicitly broken by a term oc fi 2 (f>\(f>2- 
Such term, for /i <C v induces a type-II seesaw among the scalar vevs of <pi and <p2 ■ (4>i) <C {4>j)i i 7^ j- As noted 
in |48j . such explicit breaking can be circumvented by the introduction of an additional scalar, say 03, whose vev 
generates the requi red t erm: = /i'(</> 3 ). It is exactly along those lines that we build our scalar potential. Indeed, 
provided that // in (4.2 1 is suppressed, fx' <C 1 GeV, the minimization of Vsc admits two possible hierarchical patterns 
for the vevs: W3 <C i?2,i and v 2 -C ^3,1- As we will show below, only the latter is physically viable. One may wonder 
about the naturalness of such a suppressed mass parameter \j! ' . Let us stress that very small values of \j! are actually 
technically natural. Indeed, this term, as well as the couplings h and y2 , are all terms linear in iJ 2 - By setting them to 
zero, one actually enlarges the symmetry group by an extra U(l) factor. Therefore, small values of these parameters 
are natural, in the 't Hooft sense ■ 



A. CP odd neutral scalars: A and 3 



Three CP odd neutral scalar fields arise from the spontaneous breaking of the electroweak symmetry: one pseudo- 
scalar Zl , the longitudinal polarization of the gauge boson Z, one massive pseudo-scalar A and the massless Goldstone 
mode, associated with the spontaneous breaking of the global symmetry U(1) B ^, the Majoron J. 7 

The mass eigenstates are obtained by the basis transformation 



= Rps\ J , (4.8) 





where Rps is a 3 x 3 orthogonal matrix 

/ cos(/3) sin 2 (/3)A - sin(/3) tan( 7 ) A 

R PS = sin(/3) - sin(/3) cos(/3) A cos(/3) tan( 7 ) A | , (4.9) 
\ tan( 7 )A sin(^)A 

with A = cos( 7 ) / v/l — cos( 7 ) 2 cos(/3) 2 and the mixing angles j3 and 7 are by definition 

tan(/3) = — , tan( 7 ) = — . (4.10) 

The angles f3 and 7 control the coupling of J to SM fermions. Indeed, the interaction term relevant for the Majoron 
phenomenology is 

-CDigjff 7l5fJ, (4.11) 

with 

g Jff = wi{fi) tan(/3) A (4.12) 



v 



and to/ is the fermion mass. Strong constraints apply on these couplings, stemming from star cooling processes |52j . 
In particular, the experimental upper limit on the cooling rate of white dwarfs implies: \gj e e\ ^$ 10~ 12 . Then, from 
pTO} and Ell we obtain in the limit (3 < 1 

P < 7 x KrVtan( 7 ) , (4.13) 



implying V2 < 0.2 GeV \ fv3Jv. As already stated at the beginning of this section, a hierarchical pattern of the type 
v 2 *C W3 < fi can be easily fulfilled, as the scale of V2 is directly related to the dimensional parameter /i': 

v 2 ~ -2-= ; — 2^ ' ( 4 - 14 ) 



7 The Majoron is exactly massless in our setup. Notice that J ', being the Goldstone boson of a spontaneously broken global symmetry, 
may acquire a mass through gravitational effects, as shown in 1501 151 |. However, we will not consider this possibility in the following. 
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which is suppressed compared to 1^3 as p! can be naturally set to a scale much smaller than the EWSB scale. The 
second physical pseudo-scalar, A , has mass 

M% = fi' ^ +«3 2 ) + tfo3 a „ , = , v cos(/3) _^% (4.15) 

A p 2 Vl v 2 v 3 2v 2 p ^2tan(/3) v ; 

The actual value of M40 depends on the ratio /v>2. Since the couplings of ^4° to the SM fermions are sin(/3) 
suppressed, Ma° and thus /i'/i>2 are unconstrained. 

B. Charged scalars: 



As in any two-Higgs doublet model, the charged scalar spectrum is composed of one physical field H ± and the 



eaten longitudinal degree of freedom Wf~- They are related to the interaction fields Hi 2 through the orthogonal 



transformation: 



H+\ _ /cos(/3) -sm(/?n f W+ 
H+ J \ sin (J3) cos [f3) J \ H+ 



(4.16) 



The charged scalar mass is given by 

a v 2 / tan(7) ^ , \ Af 2 W 2 , 

M ^ = y UnWV " Kl2 J - " T Kl2 ' (4 - 17) 



where in the last expression we used the approximation given in (4.15 ). Since Mfj± > 0, one requires k' 12 < 2 M^ /w 2 . 
An experimental lower bound on M H ± is obtained from iJ* pair production at LEP [53 and the subsequent decays 
H — > t v T and H ± -> cs. For to#± < mw, H ± decays only to SM fermions, so the bound M^± > 78.6 GeV 
applies. 

C. CP even neutral scalars: h° , H° and Ha 

We introduce the CP even mass eigenstates h°, H° and Ha, which are related to the interaction fields ft.1.2,3 through 
the basis rotation: 




(4.18) 

In the limit /i',«2 <C v 3 < fi, the CP even scalar mass matrix can be further simplified and Rns just consists in a 
rotation of angle 9 between the eigenstates h° and H°. Moreover, at leading order in /?, Ha and the pseudo-scalar 
A are degenerate in mass and both decouple from the other particles, so no constraints apply on Mh A ■ Within this 
approximation and introducing as a shorthand 

mi = vi y/2 Ai , m 3 = v 3 \/2X 3 and m 13 = y/v{ v 3 k 13 , 
the masses of the neutral Higgs H° and h° and the mixing angle 9 are given by the relations 



M 2 H/h = 



1 / 1 7, \ / ml — m? + \ (m\ — m?) 2 + 4mf, \ 

- (m 2 + m 2 ± yV 2 - m 2 ) 2 + 4mf 3 J , 9 = Arctan \m\ 3 ' (419) 



and Myo > M/jO. The mixing angle |0| takes values from zero to ir/2. The couplings of H° (h°) to SM particles are 
given by the SM Higgs ones times cos(#) (sin(0)). For maximal 9 ~ tt/2, h° couplings are unsuppressed compared to 
the SM case, so LEP-II bounds apply and M H o > M h o > 114.4 GeV [54 . In the opposite case, with suppressed mixing 
angle \9\ -C 1, only H° get sizable couplings to the SM, and the former bound on M H o still applies. Conversely, for 
16*1 <C 1, LEP-II bounds are rather weak in constraining the mass of the lightest Higgs. Notice that h° contributes to 
the invisible Z decay; however, the Z—hP — J coupling is /3 4 suppressed: this contribution is negligible and no relevant 
constraints apply on M h o from this decay. Nevertheless, for sin 2 (0) > 0.1, LEP-II bounds imply M h o > 80 GeV. In 
the following, we assume the conservative limit M^o > 114.4 GeV, which is valid for all values of 8. 
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An almost invisible Higgs boson 

As occurs in models with multiple scalars, the Higgs bosons may decay invisibly. In our scenario, both H° and h° 
can decay into two Majorons, thus precluding their detection at present particle colliders, LHC included. 
The total decay widths of H° and h° are given by 



T(H/h) 



-E 

8tt ^ 



K (M H/h ,M i ,M j ) 
™l/H 



\M l3 \ 2 e(M 2 H/h -{M l + M 3 f) 



(4.20) 



In the equation given above the kinematical factor is k(x, y, z) = J (x 2 — (y + z) 2 ){x 2 — (y — z) 2 ). We consider below 
for simplicity only tree-level two-body decays into identical final states. The decay probabilities of H° and h° to 
neutral scalars are proportional to the norms of the trilinear couplings, which at zeroth order in (3 read: 



M 



H/h 
J 



M 



H/h 



M 



M 



H/h 



H/h 



H/h 
J J 



H I h 

' V>7 = C0S (/3) K 13 



cos(0) 
sin(0) 



2cos(/3) tan(7) A3 



.H/h 

A H+H- 



,H/h _ 
A H+H- ~ 



cos(/3) K 



12 



cos(0) 
sin(0) 



cos(/3) tan(7) K23 



sin(6») 
-cos(6») ' 

sin(6») 



cos(0) ' 



= v \ a aa \ > A 



H/h,2 \H/h _ 



A A 



H/h 
SS 



.H/h 

' A ss 



cos(/3) ki 



cos(/3) T\ 



cos(<2) 
sin(0) 

cos(0) 
sin(0) 



cos(^) tan(7) k 23 



J3 cos(/3) tan(7) 



sin(0) 
-cos(0) ' 



sin(0) 



(4.21) 



(4.22) 



(4.23) 



(4.24) 



The decay probabilities of H° and h° into SM particles are similar to the SM ones, see e.g. [S5 
on the mixing angles 9 and /3. At tree level, for the decay probability into fermions, we have 



-cos(0) 
modulo a dependence 



M 



H/h 
f 



2 N, 



M 



H/h 



1-4 



M H/h, 



COS 2 (0) 

sin 2 (0) 



(4.25) 



where N c is the number of colors and m/ is the fermion mass. The tree-level H° /h° — > W + W decay probabilities 
depend on 



9w 2 

— — v 
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\w /h (M H /h) 



^w h ( m ) = cos(/3) 



m 



Mi 



w 



cos(0) 
sin(0) ' 



(4.26) 
(4.27) 



where M\\r is the VK-boson mass and gw is the weak gauge coupling constant. A similar expression holds for H /h 
decays into pairs of Z bosons. 

From the expressions above, we can estimate the invisible branching ratio of the Higgs bosons. First of all, obviously, 
the smaller the quartic portal couplings, the smaller the invisible Higgs decay widths. Second, the decays into Ka 
or A , even if equal at leading order in f3, do not have a similar impact on Higgs searches at colliders. Indeed, the 
pseudo-scalar A eventually decays almost exclusively into SM fermions: H° /h° — > A A can thus be considered as 
a visible channel. Conversely, Ha essentially decays into Majorons. 

In the low mass region, M H / h < 2Miy, neglecting the masses of the decay products, the Higgs invisible to visible 
decay width ratios are 



T(H/h -> inv) T(H/h -> J J) + T(H/h -> h A h A ) + T(H/h S S) 



T(H/h -> SM) 



T(H/h ^bb) + T(H/h -+ A° A°) 



(4.28) 



provided Mh A — M A a < Mjf//,/2. The decay width to 6-quarks is Yukawa suppressed, so the ratio above is simplified 
to 



T(H/h -> inv) _ T(H/h -> J J) + T(H/h -> S S) 
T(H/h ->■ SM) ~ 1+ T(H/h — >• A A Q ) 



(4.29) 



Consequently, for low masses, both H and h mostly decay invisibly. 
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H° dominant decay channel 
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FIG. 11: Dominant H° (top) and h° (bottom) Higgs decay channels in function of their masses. 



In the high mass regime, Mjj/h ^ 2Afvv, under the approximation that the Majoron channel constitutes the main 
invisible decay and the visible channel is mostly due to decays to gauge bosons, we have 



T(H/h ->■ J J) 



16 M; 



r(#/ft inv) ^ 
T{H/h -»• SM) ~ T(H/h -^W+W-)+ T(H/h -> Z Z) * ^ M 4 



Ki3 T 2 A 3 tan(7) 



tan(0) 
l/tan(6») 



(4.30) 



From the previous estimate we infer that, for a maximal 8 <~ tt/2, the heaviest Higgs boson, decays prevalently 
into two Majorons, thus forbidding its detection at current collider searches. The opposite occurs for the lightest CP 
even scalar hP . On the other hand, for higher values of Mho (M^o) and a sufficiently small (large) mixing angle 9, 
the visible decay rate of H° (ft ) becomes sizable. It dominates for very heavy Higgs bosons. 

we display the frequency at which the H and ft decays channels are the dominant ones, displayed in 



In Fig. 
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the top and bottom panels, respectively. In order to produce this plot, we use the Higgs decay branching fractions 
computed by the program micrOMEGAs [56], that we also use to study the Dark Matter sector, as discussed in 
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Section[V] As expected, we see from Fig. 11 that above the W threshold, the heavier the Higgs bosons the larger their 
visible decay rates. 8 Conversely, in the low mass regime the Higgs bosons are clearly unobservable as we explained 
above. 

V. DARK MATTER 

We discuss in this section the third building-block of our model: the existence of a viable Dark Matter candidate. 
Below the EWSB scale, the complex scalar S is split into two real components So and Si, the lightest one being the 
DM. Real scalar singlets provide the simplest DM candidates, for which a large literature exists [57]. In our model, 
we shall stress two important aspects: first, the stability of DM is not an ad-hoc prescription, but results from the 
remnant Z 2 , So or Si being the lightest particle odd under this discrete symmetry; second, we emphasize again that 
introducing the scalar S not only provides a DM candidate, but is also necessary in our leptogenesis scenario. 

The masses of the two real components of S are: 



-! afl) = Ms+o (^i«i+^2«2+^3«!) ± (hv lV2 - n"v 3 ) . (5.1) 



(l) r~o 2 



The mass splitting in this case is controlled by the parameters h and fj," . However, since v 2 <C W3, the latter term 
dominates and mgo < mg 1 for positive /j," . 9 



As seen from Vdm, eq. (4.3), S has several portal couplings to the Higgs fields, implying many annihilation chan- 
nels [SSj. Like in most of the singlet scalar DM scenarios, S easily gets a thermal relic abundance in agreement with 
cosmological requirements. 

A. Relic density 

The DM annihilation cross-section can generically be written as 

A 2 

(5.2) 

m S 2 

where the effective coupling A e // indicates that each annihilation channel receives in general several contributions. 
When S annihilates into scalars, the cross-section is the (coherent) sum of the contact term interaction, for which 



\ e ff J-'i, cf. eq. (4.3), and of scalar-mediated interactions, where \ e ff depends on the different trilinear scalar 

couplings, such as \ s s introduced in the previous section. 

For light DM, that is 7715 < Mw, S mostly annihilates to Majorons, as well as to pairs of Ka or A , granted the 
latter are light enough. Notice that the annihilation cross-sections into pairs of Ha and A coincide at zeroth order 
in (3. 

For heavier DM, new annihilation channels are open. In the case > Mw, S can annihilate into pairs of 
through h° and H° s-channel (the Ha mediation is (3 2 suppressed): 10 



9w ( VAf s A&(2m s ) X h s s \ h w {2m s ) 



16 7r m| V Mjjg - 4 m| M% - 4 m| 



(5.3) 



where \^J h {m) were introduced in eq. (4.27). A similar expression holds for the annihilation into pairs of Z bosons. 



In the high mass range, S may also annihilate into pairs of charged H ± , or to pairs of CP even scalars h° and H°. 



Increasi ng th e DM mass, the quartic couplings Ti which control the DM mass, eq. (5.1 1, and the effective couplings 



Af s \ eq. (I424J), should increase as well so that the annihilation cross-section remains large enough, in order to obtain 



the observed DM relic abundance. 



8 Notice that we only consider the two-body decay widths H°/h° — > W + W~ . However, in the SM the tree-body decays through off-shell 
W actually dominate for M H j h > 135 GeV, cf. | 55| . in which case the Higgs visible decay channels should prevail here as well. 

9 In the following we will however denote by S the DM candidate. The hea vier state will decay to DM plus Majoron. 

10 For simplicity, the widths of the Higgs fields have been neglected in eq. | |5.3| , although they are taken into account in the numerical 
evaluation. 
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Dominant DM Annihilation Channels 
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FIG. 12: Main Dark Matter annihilation channels for different Dark Matter mass ranges. 



Numerical evaluation 

In order to accurately determine the relic abundance of S, we implemented our model in micrOMEGAs [55] , through 
the program FeynRules |40j . We then performed a scan over the full scalar parameter-space, by assigning random 
values to the different couplings. All A and k quartic couplings were varied from 10~ 4 up to the perturbative bound 
47r. The trilinear coupling h was chosen between 10 -6 and 10~ 2 . The scalar masses were randomly varied from their 
experimental lower bounds, discussed in the previous section, up to 500 GeV. In particular, as regards the CP even 
scalar masses, recall that we impose the conservative bound Mjjo/ho > 115 GeV. We vary the mixing angle 9 in the 
range: < \8\ < ir/2. For the unconstrained scalars Ha and A , their (almost degenerate) mass was varied between 
1 GeV and 100 GeV. 

The trilinear mass term p!' was scanned over in the range (1-10 2 ) GeV, while // typically took values between 10 
eV and 10 MeV. Finally, ps was varied from 1 GeV to 500 GeV. 

We demand the relic density of S to account for all the DM abundance and to lie within the 3cr range of WMAP 



.10: 



DM 



0.229 ±0.045. 



We illustrate the relative contributions of the different annihilation channels in Fig. [12] Binning the DM mass range 
into intervals of interest, we present the frequency at which a given channel is the dominant one. For example, before 
the W channel is kinematically open, i.e. for ms < Myy, we see from Fig. |12| that S annihilates only into pairs of 
J , A and Ha- For heavier DM mass, new annihilation processes are possible. In particular, annihilation into gauge 
bosons, charged scalars or CP even scalars tend to be the dominant processes. Notice that Fig. [12] only displays the 
frequency a given annihilation channel dominates in a given mass interval and not the relative contributions of the 
different channels. 

B. Direct detection constraints 

The Dark Matter can scatter on nucleons through scalar-mediated i-channels: the spin-independent (SI) elastic 
cross-section receives contributions from both h° and H exchange, according to: 



1 /4,n 

47T m 2 So 



2 j-2 
m nfn 



X 



H« 



Mho 



Ml, 



(5.4) 
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FIG. 13: Spin- independent cross-section against ms- the blue points are the model predictions which provide the required relic 
density. The red line represents XENON100 results, extracted from [21] , 



In this expression, Hs,n is the S'-nucleon reduced-mass and m n the nucleon mass. The factor /„ is the effective 
Higgs-nucleon interaction and varies from 0.14 to 0.66 [59j . The couplings A#o and A^o are given by 

\ - 1 x H/h / cos(6») , . 

XH/h -^(f) Xss \sm(6) (5 - 5) 

Assuming the conservative bound M H o/ h a > 115 GeV, we see from the previous expression that for 9 w (tt/2) the 
main contribution comes from H° (h°) exchange and a^ 1 is then mostly affected by JF\. Notice that, contrary to |60) . 
where the mixing suppression 9 <C 1 was balanced by a very light scalar h° (M^o < 1 GeV), in the present scenario, 
taking M h o above the LEP-II bound drastically forbids such an enhancement. In the limit of small mixing angle 9, 
assuming /„ ~ 1/3, the SI elastic cross-section can be roughly expressed as 

<■ ^ s fl(-^)~s* ~ (^y' (g) 2 , (,e, 

which shows that S can easily saturate current direct-detection bound for electroweak scale DM [T5]-[2T]. As we saw 
in the previous subsection, since the annihilation cross-section scales as (Tk/ms) 2 , the couplings T% should be sizable 
for large DM masses, otherwise S relic density would overdose the Universe. This, in turn, implies that for heavy 



DM the scattering cross-section on nucleons is almost independent of the DM mass, cf. eq. (5.5 1. 

The dependence of a^ 1 on mg in the low and high DM mass regimes is manifest in Fig. 13 In this plot we compare 
the model predictions (blue dots) for <r„ with XENON100 results [3T] (red curve). We can see that while only a small 
region of the parameter-space is already excluded by current data, the next generation of direct-detection experiments 
would probe a large part of it [61j . Notice in particular that in the low mass regime high cross-sections can be reached, 
due to non-suppressed J^i^ couplings. A light DM with large J-1,3 couplings is possible through a partial cancellation 



in eq. (5.1 1, which depends on the value of the parameter fi". 



VI. CONCLUSIONS 



In this paper we study a seesaw extension of the Standard Model based on a global U(1) B _^ symmetry group, where 

L can be thought as a generalized lepton charge. This global symmetry is spontaneously broken at the electroweak 
scale. Suitable scalar and fermion representations are added to the SM particle content so that a tiny Majorana mass 
for active neutrinos is naturally generated, in agreement with neutrino oscillation experiments. More specifically, an 



21 



extra Higgs doublet H 2 and a Higgs singlet H 3 are added to the SM, together with a heavy Dirac fermion Nr>. The 
lepton doublets and Njj interact through neutrino Yukawa couplings which can violate the lepton number. When Np 
mass is set at the TeV scale, the model realizes a UV-complction of the inverse-seesaw mechanism. 

We show that, with the addition of two extra SM-singlets in the model, a Majorana fermion N 3 and a complex 
scalar S, it is possible to explain quantitatively both the observed baryon asymmetry of the Universe through an 
original leptogenesis mechanism and the Dark Matter relic abundance. 

Leptogenesis in this model is implemented in two steps: first an asymmetry in N D and S is generated by the 
out-of-equilibrium decays of N 3 . In a second step this asymmetry is converted in a non-zero lepton charge due to 
fast neutrino Yukawa interactions. The latter constitute a link between leptogenesis and neutrino mass generation. 
We solve numerically the Boltzmann equations relevant for this two-step leptogenesis scenario and show that the 
observed amount of baryon asymmetry is easily achieved. We concentrate the discussion on a TeV scale scenario, 
and show that, provided neutrino mass constraints arc fulfilled, no lower-bound on the mass of Njj is imposed by the 
requirement of a successful leptogenesis. However, this scenario of leptogenesis is also viable at much larger scales. 
An important feature of this mechanism is that the source and damping terms do not depend on the same couplings, 
therefore large CP asymmetries can be obtained even in the regime of weak washouts. 

In the second part of the paper, we analyze in detail the mass spectrum of the model and provide constraints on 
the parameter-space arising from low-energy physics. In particular, we show that the presence of a massless Majoron, 
which corresponds to the Goldstone boson associated with the spontaneous breaking of the global U(l) B _i symmetry, 
has an important impact on Higgs boson searches. Indeed, light Higgs scalars H° and h°, M H o/ h o < 140 GeV, would 
mainly decay into pairs of Majorons, thus making difficult their observation at colliders, LHC included. 

Concerning the scalar field S, after the breaking of the electroweak symmetry, the lightest component of £ remains 
stable, due to the presence of a remnant Z 2 symmetry, and provides a viable candidate for Dark Matter. Its mass is 
unconstrained and can take values as light as few GeV up to few TeV. Numerous annihilation channels are present, 
allowing the relic DM density to be consistent with cosmological observations. We study the possible signatures 
of DM in direct detection experiments, and show that while the current constraints exclude already a part of the 
parameter-space, a large region may be probed by the next generation of detectors. 

This model explains in a common framework three main experimental issues: neutrino mass generation, the baryon 
asymmetry of the Universe and the Dark Matter relic density. Many observables are predicted, but their measurements 
probe uncorrelated sectors, making this minimal extension difficult to falsify. 
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Appendix A: Interaction Lagrangian of the Model 

The full interaction Lagrangian of the model is: 

£ int - £f™ - V SB - V DM - MN D N D - ^M 3 N 3 N% (Al) 
- (y[N D Hl£i + y 3 2 N c D HU 3 + ^H 3 N D N C D +gSN D N 3 + h.c.) , 

where Cf^ is the Yukawa interaction Lagrangian of the Standard Model and 

Vsb = -Mi Hi H x + X 1 (Hi H x ) 2 - fj 2 Hi H 2 + A 2 [Hi H 2 f - fx 2 H* 3 H 3 + A 3 {H* 3 H 3 ) 2 
+ k 12 Hi H X H\ H 2 + k' 12 H\ H 2 h\ H x + k 13 H\ HtH* 3 H 3 + k 23 h\ H 2 H^H 3 

- -]j= (Hi H 2 H 3 + Hi HiHfj , (A2) 
V DM = n 2 s S*S + \ s (S*S) 2 + ^1^1^15*5 + T 2 HlH 2 S*S + T 3 H* 3 H 3 S*S 

+ h S 2 Hl H 2 + h* S* 2 H\ Hi - ^-{S 2 H* 3 + S* 2 H 3 ) . (A3) 

v2 
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Appendix B: Computation of the CP Asymmetry 

The relevant interaction Lagrangian which is involved in the generation of the CP asymmetry in the out-of-equilibrium 
decays of the Majorana neutrino N 3 is the following: 

-C mt D ^=H 3 N D N D + gSN D N 3 - ^S 2 H* + h.c. , 

— T — T 

where N D = CN D , N 3 = N 3 = CN 3 . The CP asymmetry in the decays of iV 3 is defined as: 



(-CP 



r (n 3 -> n d + s) - r (n 3 -» n d + s) 

T (N 3 -> N D + S) + T (7V 3 -> N D + S) 



(Bl) 



where E{n} indicates the sum over all possible on-shell states in the loop of Fig. jlj while the phase space factor in 
the integral is, in general 



dIL = d3pni 



d 3 Pn k 



(2n)32E ni (2nf2E ni 



(27r)V 4 ) [p N3 - > k ^ 2 - 



(B2) 



Pat 3 and p Uj (j = 1, ...,k) being the 4- momenta of the decaying Majorana neutrino N 3 and the final state rij, 
respectively. We consider the physical intermediate processes: 11 7V 3 — > Np + S and Njj + S —> Njj + S. The 
corresponding tree-level amplitudes read: 12 



iM (0) (N 3 -> JV D + 5) = Z. 9 ^(pAr)C , - 1 UAr3(p A r 3 ) : 
iM (0) (N 3 -> 7V D + 5) = -^U A r(p^)M Ar3 (p A r 3 ), 



(B3) 



iM (0) (7V D + N D + S) = -- 



ft 



Ph 3 



-mi 2 



(a*v^(p N )C 1 u N (p' N )) 



We perform the product of the three amplitudes in (B3) according to eq. (Bl I and sum over the polarizations of the 
outgoing fermions. After some algebra, we get 

M {0) {N 3 ^N D + SyM m (N 3 -> N D + S)M {0) {N D + S^N D + S) = 



2 ^ n"M 2 M 3 \ i ( PN ■ P ' N ) ((p N +j/ N )-p Na ) 



2g a 2 



M 2 



MM 3 



Integration over the phase space 



The relevant integrals in the numerator of |BlJ) are 

d 3 p' N d 3 p' s S„ 



In 



(2^2E' N {2irf2E> s p% 3 
where S n £ {{p N -p' N ), (p' N -Pn 3 ), (pn -Pn 3 ), M 2 }. 



(2n)^(p N3 -p' N -E' s ) 



(B4) 



11 Notice that the other possible cuts in Fig. ^do not contribute to the CP asymmetry as they do not correspond to physical processes. 

12 In the following we indicate with 7713 the thermal mass of the scalar singlet H3, which provides an infrared regulator of the the N3 decay 
one-loop diagram. 
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It is convenient to express /„ in terms of adimensional quantities, mainly: a = En/M 3 , b = \pn\ /M 3 , x = M/M 3 , 
xs = m s /M 3 and x 3 = m 3 /M 3 : 



h 



^2,3,4 



2x 2 — x 2 

2 k(1, x, x s ) H t — - C (x, xs, x 3 , a, b) 



1 

32^ 

{B{x,xs),2a, 2x 2 } C (x, x s , x 3 , a, b) 



(B5) 
(B6) 



where k is a kinematic factor introduced below eq. (4.20 ), B(s, t) = y k (1, s, t) + 4s 2 and 

'2s 2 -u 2 -aB(s,t) +6/^(l,^ ) s) , 



C (s, t, u, a, b) = log 



2s 2 - u 2 - a B (s, t) - b k (1, t, s) 



(B7) 



Now we complete the integration over the phase space in the numerator of eq. |B1J). The relevant integrals can be 
arranged in the form: 



The full computation results in: 



1 

4?r 



da — /„ S J 1 — a — y a 2 — x 2 + x 2 s 



(B8) 



Ji 



J-?, — J 3 — 



J A = 



1 



1287T 2 

1 



-k(1,x,xs) + (2x 2 — x 2 )\o§ 



B (x, x s ) log 



x\ + K (1, X, XsY 



128^ 2 y '" , '" 0/ '" \x 2 + K (l : x,x s ) 2 



x log 



64tt 2 



x\ + K (1, X, XsY 



A similar computation applies for the denominator of eq. (Bl). We have in this case: 

2 



dU 



N.S 



M 2 

M {0 HN 3 ^NS) = g 2 -^K(l,x,x s ){2x + B(x,x s )} 

47T 



(B9) 
(BIO) 
(Bll) 

(B12) 



The CP asymmetry in the decays 



Taking into account the results obtained in eqs. (B9)-(B12) and the general expression (Bl), we get the final 
expression of the CP asymmetry: 



1 ^" t 1 \ F 2(x,x s ,x 3 ) + 2x [x + B(x,x s )]Fi(x,x s ,x 3 ) 
ecP = "16^M3 (a) 2x + B(x,x s ) 



(B13) 



where 



log 



F\ (x, xs, x 3 ) — . , , 1 . 

k{1,x,x S ) \xj + n(l,x,x s ) J 

F 2 {x,x Sl x 3 ) = -k (l,x, x s ) + {2x 2 - x 2 3 )Fi{x,x Sl x 3 ) . 



(B14) 
(B15) 



Therefore in the limit mj, M <C M 3 , which we are interested in, we get the approximation reported in eq. (3.2) 
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Appendix C: Boltzmann Equations 



In this appendix, we introduce the set of Boltzmann Equations (BE) that are used for the numerical evaluation 
of the baryon asymmetry. More details on the network of BE can be found in the appendices of [35] and [44] . For 
a given particle (asymmetry) X, we denote as usual by Yx its comoving number density, i.e. the number density 
normalized to the entropy density. We assume Maxwell-Boltzmann statistics for both fermions and scalars. In an 
expanding Universe, the evolution of Yx is governed by the Boltzmann equation: 

dY 

3H{z) ^t =_ ? [Xav±ij] , 

a,i,j, ... 



where 

Yx Y Y Y 

[XaV±ij] = -^-^j eq (Xa^ ij) - -k-^j eq (ij -> X a) 

1 X la X i 1 n 



and z = M 3 /T is the evolution parameter, while 7 e? are the equilibrium reaction densities of the different processes. 
We will limit our analysis to 1 -O- 2 and 2h2 processes, but will include the on-shell part of some 2 -o- 3 scatterings for 
consistency. If these processes conserve CP, then we use the notation [X a <H> ij], as -f eq (X a — > ij) = j eq {ij —> X a). 
In a radiation dominated Universe, the Hubble constant H(T) and the entropy density s are given by 



An* g, T 2 2tt 2 3 

H{T) = V^5-^' S = 9 *^5 T • 

In these equations, 5* is the number of relativistic degrees of freedom present in the thermal bath at the leptogenesis 
time scale. In the case of the SM, at temperatures above the EWSB, one has g^ M — 106.75. Assuming that 
the non-SM scalars S, H 2 and H 3 and the Dirac fermion Njj are relativistic particles at T w M 3 , we obtain: 
.9* = gf M + 46/4 = 118.25. 



As already explained in Section III the main source of Njy and S asymmetry production during the first stage of 
leptogenesis are the CP- violating decays and inverse decays of A3, 



l eq (N 3 ^N D S)^ lD f 1 ^^) - l eq (N z ^N D S) , 



where €qp is the CP asymmetry in the decays, defined in eq. (Bl |, and jp> is the CP conserving total decay width 
of A 3 . The last equality results from CPT invariance. 

We further include in the BE AAq = AS = 2 scatterings shown in Fig. [2j whose corresponding collision rates are 
denoted as: 

^i(N D S^ N D S) = rA2 and ^ q (N D N D O S S) = 7 1 2 . 

Note that, as in standard leptogenesis, Njj S <=^ Njj S processes mediated by A 3 in a s-channel develop an on-shell 
part, which is CP-violating. To avoid double-counting of this resonant part, already accounted for by the inverse 
decays, the on-shell contribution should be subtracted from the full Ad S «-» JV^ S scattering rate. 

In addition to the standard source term given by the decays of A3, we include the CP violation arising from the 



2f>2 scatterings involving an external A 3 , which also depends on the CP-violating phase a entering in tap, eq. (3.2 1. 
The CP asymmetry for each diagram is computed as in the standard leptogenesis scenario, e.g. [B2] and [63 . However, 
in our model a contribution to CP asymmetry in the A3-scatterings arises from both s-, t- and it-channels, as depicted 
in Fig. [3] The corresponding thermal rates, in this case, are: 



a) 7 e9 (A D A 3 <± H 3 S) = 7^3 (1 ± £cp) = l eq {N D N 3 ^H 3 S). 

b) 7 e "(A 3 S^N D H 3 ) ee 7^ (1 T e b CP ) = 7 e9 (^3 S t* N D H 3 ) . 

c) ^{Nd S v± A 3 H 3 ) = 1 C N , (1 T e c C p) = l eq (N D S <=> A 3 H 3 ) . 

d) 7 (A 3 S^Hxi) ee 7 ^ 3 (1 ± e CP ) = 7 (A 3 S^^l). 
c) 7 (A 3 S<=tH 2 l) ee 7^ (1 ± e C p) = j(N 3 S U H 2 £) . 
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The CP asymmetries in the scattering, t k CP (k = a, e), are defined by e^p = ecp AK k , with 
AK a, b = ( 7 ^) t -(^3) 8 (7^) t -(7^) u A d = AK e = 1 

a,b ' _,c ' 

7jV 3 7at 3 

Here (7/v 3 ) c j c = (s,t, u), corresponds to the s-, and u-channels of the different processes shown in Fig. [3J Notice 
that, similarly to the ANo — AS — 2 scatterings considered before, as explained in [S3], we have to subtract the 
resonant CP-violating contribution of the 2f)3 processes in which N 3 is exchanged in s-channel. The non-resonant 
parts of such processes are not taken into account in our computation, since they are at higher order in the couplings. 

As regards N 3 three-body decays, which are at the same order in the couplings as2f>2 scatterings on N 3 , they 
are phase-space suppressed and so give a sub- leading contribution with respect to the two-body decays |63j . and we 
consequently do not include them. 

We further consider the effect of S self-annihilations (see Fig. [4 1 , which could washout the asymmetry Y&g for large 



values of the coupling h (see eq. (4.3 1). The related interaction density rate is noted 

7 eq (SS o H 1 H 2 )= 1SS . 

Several processes participate in the second phase of leptogenesis. Besides the scatterings on N 3 , the 7^ discussed 
above, we include the following interactions, at the lowest order in the neutrino Yukawa couplings: 

a) No decays and inverse decays: "f eq (N D <-» l a H\) = 'foe an d ^ eq {N D «-> £pH 2 ) = l D j- 

b) AL = 1, iJi-mediated scatterings with top-quark: j eq (No£ ■<->• Q 3 i) = 7jy D for the s-channel contribution and 
"f eq (Noq3 O It) + ^ eq (N D t o lefe) = 27^ for the t-channel contributions. 



As already stated in Section III the leptons participate in No mediated AL = 2 scatterings: 7^ = ^ eq (th\\ -H> £H 2 ) 
and^Ef H 2 ). 

We are ready now to report the complete set of Boltzmann equations relevant for the computation of the baryon 
asymmetry of the Universe in the two-step leptogenesis scenario described in the text. We include all the in- 
teraction terms introduced above and we use the simplified notation y^ 3 = YN 3 /Y^ q , y&x = Yax/Y^ 1 and 
Y' x = (s H{z)) dYx/dz. At first order in the asymmetry (zeroth order for N 3 ), the full system of Boltzmann equations 
is the following: 

Y k = - l N 3 N D S] - [N 3 ^±N D S]- [N D N 3 <± H 3 S\ - [~N D N 3 <± H 3 S] 

- [SN 3 f=i N D H 3 ] - [SN 3 «=* N D H 3 ] + [S N D <=i N 3 H 3 ] + [SN D <± N 3 H 3 ] 
+ [£Hi <=* N 3 S] + [lUi r± N 3 S] + [H 2 £^N 3 S] + [H 2 1t±N 3 S] , 

Y' Nd = [N 3 t±N d S]-[N d Sv±N d S]-[N d N d ^ SS] 

- [N D N 3 +i H 3 S] - [N D H 3 t± N 3 S] - [N D S^N 3 H 3 ] 

- [N D o tHj] - [N D o lH 2 ] - [N D l^q 3 i] - [N D q 3 o It] - [N D t^l7E\ , 

Y s = [N 3 ^ SN D ] - [SN D N D ] - [S S o N D N D ] 

- [SH 3 ^ N 3 N D ] - [S N 3 «=i N D H 3 ] - [S N D A^ 3 H 3 ] 

- [SN 3 <± £Hx] - [SN 3 *=t £H 2 ] - [SS^ HrH 2 ] , 

Y't = [No O tHt] + [No O tH 2 ] - [I N D Hftt]-[it« N D q 3 ] - [i© O No t] 

- [eH^NaS] - [£H 2 ^N 3 S] - [tHi**TH 2 ] - [££^H 1 H 2 ] , 

Y' Hi = [N D ^H 1 i\-[H 1 i^N 3 S}- [H l H 2 ^ S S] - [IH X ^IH 2 ] - [H l H 2 ^II\ , 
Yh 2 = Wo O IH 2 ] - [H 2 l z± N 3 S\ - [H 2 Hi «■ SS] - [IH 2 o IHx) - [H x H 2 o 11] , 



2G 



V' 

Y H 3 



[H 3 S^N 3 N D ] - [H 3 N D «=* N 3 S] - [H 3 N 3 T± N D S] 



The evolution equations of the antiparticles are obtained by taking the CP conjugates of the different rates. The 
Boltzmann equations of the number density (asymmetry) finally read: 



V' 

N 3 



(1 - VN 3 ) ID 



(CI) 



V' 

Y AN L 



(Vn 3 - 1) (ecp ID + 2 e a CP l% a - 2 e b CP 7^ - 2 e c CP 7^) 

2 (7A 2 + 27A2) (VAN D - VAS) ~ 7jV 3 (VNsVANo ~ VAH 3 " 
l b N 3 {VAN D ~ yAff 3 + VN 3 VAS) ~ 7W 3 (2/AJV d ~ VN 3 VAH S ^ 
7D-« (l/AW D - VAl - VAHt ) - 1 D J {VAN D + VAl + VAH 2 ) 
{lN D + 2 7at 3 ) iVAN D ~ VAt) , 



- VAS) 
VAS) 



(C2) 



= - (VN 3 - 1) (ecp + 2 e a CP 7^ - 2 e b CP 7^ - 2 e^ P 7^ + 2 e CP 7^ + 2 e CP 7^ ) 
- 2 (7a 2 + 2 7A2) (z/as ~ Uan d ) - 2 7ss (2y A s - J/AHj + Vah 2 ) 



1n 3 (vas - Vah 3 + VN 3 yAN D ) - 7w 3 {vn 3 Vas 



1n 3 {VAS - VN 3 VAH 3 + VAN L 
In, (VN 3 yAS + UAH 2 + VAi) 



7jv 3 (vn 3 Vas 



VAH 3 
VAHi 



VAN D ) 
VAi) 



(C3) 



v 

1 At 



^Aifi 



= - (2/jv 3 - 1) ( 2 £ cp 7at 3 - 2 £ cp 7at 3 ) - 7d 1 {VAi + Vah x - Van d ) 

- 1 D I (VAi + VAH 2 + VANo) + (l S N D + 2 lN 3 ) (VAN d ~ VAl) ~ (lu + 2 7«) ( 2 2/A£ + VAHi ' 

- 7w 3 (2/AJ?! + 2/A£ - Vn 3 Vas) - 1n 3 (vah 2 + Vai + Vn 3 Vas) , 

= ~ 2 (?/W 3 - 1) e CP 7iV 3 - ID i {VAHi + VAl - VANo ) ~ ("fa + lu) (2yA£ + VAHt + VAH 2 ) 

~ lN 3 (VAH 1 + VAl - VN 3 VAS) - 1SS (VAHt ~ VAH 2 ~ 2 2/As) , 

= 2 (VN 3 ~ 1) e C P 7AT 3 - iDl (VAH 2 + VAl + VANo) ~ (lit + 7«) ( 2 yA^ + VAH 1 + VAH 2 ) 

~ 7 Ns (VAH 2 + VAl + VAS VN 3 )lSS (VAH 2 - VAHi + 2 VAs) , 

= (VN 3 - 1) ( 2 e cp Yk ~ 2 £cp 7^/ 3 - 2 £cp 7at 3 ) - 7w 3 (vah 3 - Vas - VANnVm) 

- 7W 3 (2/AH 3 - VANo - VASVN 3 ) - 1n 3 {VAH 3 VN 3 ~ VAN D - VAs) ■ 



VAH 2 



(C4) 
(C5) 
(C6) 
(C7) 



Appendix D: Chemical Equilibrium Conditions 



M-Hi 

US 
fJ,N D 



CB-L 



J_ 

16 



11 
16 



16 



CjVr 



15 
Hi 
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16 
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16 



CS 
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cb-l 



16 



11 
16 



21 
' 16 



16 



CjV r 



1H 
i« 



101 
16 



231 
16 



16 



Cs-L 



14 



14 



14 



TABLE II: Chemical equilibrium coefficients in eq. (D12 l for the three cases discussed in the text 



We derive in this section the chemical equilibrium conditions provided by all the interactions which are in equilibrium 
at the leptogenesis epoch, T ~ M 3 < 10 5 ^ 6 GeV. 
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The chemical potentials of each generation of SU(2)w quark doublets, Qi, and singlets, um and dju, are denoted 
by /j,Q i = (Aq, Hum = M« an d Hd m = Hd, respectively. Concerning the lepton fields, we define for each flavor a 
the corresponding chemical potentials as: Hl a = Hli He Ra = He- We denote with hn the chemical potential of Njj. 
Analogously, for each scalar field in the model we define, in a consistent notation: HHt 2 3 and hs- We remark that 
the chemical potentials of the SM fermions are assumed to be independent of the generation index, because of the 
the rapid flavor mixing interactions which occur at the leptogenesis time [42j . 

The number density asymmetries are related to the particle chemical potentials through the relations: 

Yax — gX n T Hx for bosons , (Dl) 
3 s 

Y^x — — Hx for fermions , (D2) 

where gx is the number of internal degrees of freedom of the particle X. The total baryon and lepton number 
asymmetries can be expressed in terms of the fermion chemical potentials: 

T 3 T 3 / 2 \ 

Yab = {2hq + Hu + Md) > y al = ( 2a*£ + Me + ^n d j ■ (D3) 

Taking into account the definitions given above, we have the following relations 42J: 

1. QCD and SU(2)w sphaleron interactions: 

2MQ - Mu - Hd = , (D4) 
3/iQ + He = . (D5) 

2. Hypercharge neutrality: 

3 {hq + 2[i u - Hd - He - He) + 2 (HHt + Hh 2 ) = . (D6) 

3. Charged lepton Yukawa interactions: 

He - Hh x - Me = , (D7) 
HQ + HHt - Hu = , (D8) 
HQ - HHt - Hd = 0- (D9) 

4. Lepton number conserving Dirac neutrino Yukawa interactions: 

Hn d - HHt - He = ■ (D10) 

5. (B — L) conservation: 

3 (2hq + Hu + Hd) - 3 (2ne + He) - 2hn d -2(hs + 2hh 3 - ^Hh 2 ) = . (Dll) 

We notice that the QCD sphaleron condition is redundant in this case, as all quark Yukawa interactions are in 
equilibrium. 

The different chemical equilibrium conditions enforce relations among the chemical potentials, which then can be 
expressed in terms of a subset of them. We set 

Hx = cb-l Hb-l + cn d Hn d + csHs, (D12) 

where we define hb-l through the relation: Ya(b-l) = Yab — Yal = Hb-lT 3 /(2s). We then distinguish three 
possible scenarios: 

A. Lepton number violating neutrino Yukawa interactions and S self-annihilation are decoupled, which corresponds 
to the set of equilibrium conditions 1 — 5 listed above. The different chemical potentials can be expressed in 
terms of the set (hb-l, Hn d , Hs)- 
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B. S self annihilations are always in equilibrium, but lepton number violating Yukawa interactions are still decou- 
pled. An additional equilibrium condition is enforced: 



2a*s - M-ffi + l l H 2 = . 
Only two chemical potentials are independent, that we choose to be (/j,b-l, I^n d )- 



(D13) 



C. All interactions listed above, as well as lepton number violating Yukawa interactions, are in thermal equilibrium 
during the leptogenesis era: 

Mw D + M-ff 2 + W = ■ (D14) 
In this case, all chemical potentials are proportional and can be expressed for example in terms of [Ib-l- 



The coefficients cx in eq. (D12), corresponding to the three cases listed above are reported in Tab. |D| In the first two 



cases the final baryon asymmetry is given by 



Yab = j Ya(b-l) - g Yan d 



(D15) 

In the last scenario, which corresponds to the case discussed in Section pTI| where all the interactions listed above are 
in thermal equilibrium during the generation of the BAU, we have: 



Yab = ~ Y A (B-L) 



(D16) 



Notice that expressions (D15) and (D16) should be considered valid up to the decoupling of Nd, i.e. for Tn d ^> H. 
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